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ON 

Spin correlation functions (up to the 3-site one) of disordered Ising 
model with the nearest neighbour interaction are calculated and in- 
vestigated within a two-site cluster approximation for both quenched 
and annealed cases. The approach yields the exact results for the one- 
dimensional system. The long-range interaction is taken into account 
in the mean field approximation. 
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1. Introduction 
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During recent years the statistical theory of disordered spin models is an 
object of great interest due to their various properties and applications. In 
this paper we will consider an alloy of Ising systems (AIS) that can describe 
Q\ | the processes of (magnetic or ferroelectric) ordering in alloys of magnets, 

magnets with nonmagnetic admixtures, solid solutions of ferroelectrics etc. 
This model is considered usually in two limiting cases, namely the case 
of annealed system (equilibrium type of disorder) and the case of quenched 
system (nonequilibrium disorder). There is no common opinion on the ques- 
tion, what type of disorder is realized in partially deuterrated ferroelectrics, 
e.g. C s(D x Hi- x )2PC , 4 1 , K(D X H\- X )2P O4. This situation requires both an- 
nealed and quenched cases to be theoretically studied and compared with 
experimental data. 

The exact results for both limits are known at the one-dimensional case 
[ |l[. Effect of dilution (noninteracting admixture) on the two- and three- 
dimensional Ising model (quenched case) have been studied by Monte-Carlo 
method [|, |. AIS can also be used as a test of approximative methods of 
statistical mechanics. It is known that the mean field approximation (MFA) 
yields quite satisfactory results for the Ising model. However it is not able 
to reproduce some essential properties of AIS [ |]] (such as percolation in 
diluted system and some differences between the properties of annealed and 
quenched systems). 

In spite of great attention to AIS and great advance in its investigation 
(see review in [g, ^|) the problem of calculation of AIS 's correlation functions 
(CFs) for the case D > 2 (within some acceptable approximation) is not 
solved up to now. In this paper we suggest the solution of this problem 
for both quenched and annealed systems generalizing the approach of [ 0] 
on disordered model. The short-range interactions are taken into account 
within the two-site cluster approximation (TC A) . In the appendix the long- 
range interactions are also taken into consideration within MFA. 
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2. Formulation of the problem 

We shall consider an Ising mixture with the Hamiltonian H({S, X}) 

-FH({S, X}) = H ({S, X}) = H x ({S})+H ({X}); (3 = l/k B T , (2.1) 
where 



H x ({S}) = *iSi + \ E K H S i S i (2-2) 



describes a quasispin subsystem on the sites of a simple lattice = 
1 • • • N) with the pair exchange interaction Ky in the site dependent field 
Kj. The set of spin variables {Si} represents a state of the spin subsystem 
(Si = ±l). 

The other part of the Hamiltonian 

#({*}) = E* + sE^ ( 2 - 3 ) 

i ij 

includes " nonexchange" pair interaction Vij and field Each site contains 
a spin of a certain sort and the sort configuration is described by the set of 
variables {Xi a } (a = 1 ■ ■ ■ f2, f2 is a number of sorts): Xi a = 1 if the site % 
is occupied by the spin of sort a, otherwise Xi a = 0. Interactions Kj, Kij, 
Hi, Vij depend on the sort configuration 

Hi = ^ ' HiaXi a ] Vij = ^ ' Vi a jfiXi a Xji3 ; 

a a/3 

= E K *a^Q!i Kij = E KjaJpXigXjp . (2-4) 

« a/3 

It should be noted that parameters of the Hamiltonian contain inverse tem- 
perature (3. We suppose the pair interactions to be short-range 



Via,j/3 = f3V a/ 3iTij; (3K ia jp = K a/ 3iTij; tt^ 



1, if i £ iij 



0, in opposite case ' 

where 7Tj denotes the set of the nearest neighbours of the site i (the first 
coordination sphere). In an appendix we will also take into account the 
long-range interactions in MFA. 

The model is considered at the cases of annealed system and of quenched 
one. At the first case an equilibrium sort configuration is realized and the 
system is described by the density matrix 

P ({S,X}) = Z- 1 expH({S,X}); Z = Sp{ S X y exp H ({S, X}) (2.6) 

as well as by the generating function 

F({K,fi})=lnZ = -pg, (2.7) 
where Q is the grand thermodynamic potential of the system. 
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The generating function allows to calculate the correlation functions 
(CFs) of the system 

(Si iai ■ • • Si n a n Xj l p l ■ ■ ■ ^j i/ 3 i ) c = 

S $ 5 5 , , . , . 

= j j j 7 H{K,lA) > (2-8) 

where Si a = SiXi a , the superscript c means cumulant averaging and 

((•••)) = Sp {s , JC} [/>({£ *})(■••)] • (2-9) 
The chemical potentials [ii a have to be found from the equations 

(Xia) = c a , (2.10) 

where c a is a concentration of the sort a spins. 

At the case of quenched type of disorder the sort configuration is fixed 
and independent of temperature, therefore thermodynamic averaging im- 
plies only a trace over spin degrees of freedom 

((■■■))h = S P{s} [ Px ({S})(---)}; Px ({S})=Z- 1 e W H x ({S}) ; (2.11) 

Z a; = Sp {s} exp( J H x ({5})). (2.12) 

In order to obtain observable quantities one must perform also an averaging 
over sort configurations 

((■■■)h = Sp {x} [p({X})(...)] , (2.13) 

where the distribution p({X}) is determined by the conditions of system's 
freezing. Our approximation is sensitive to the following moments of this 
distribution: 

{X ia Xjp) x = w af3 (j £ 7Ti); (X ia ) x = c a = w af 3 . (2.14) 



Free energy of a quenched system is defined as follows 

F=-p- 1 (F x ) x -0- 1 (H({X})) x ; (F x ) x = (In Sp {5} exp^ ({5})]> x . 

(2.15) 

Function {F x ) x = (F x ) x ({«}) that is the cumulant generating function in 
this case: 

m i 1 a 1 ---i n a n = ((Siiai ' " " Si n a n ) C H ) x = — — (F x ) x ({«}) • 

(2.16) 

3. Quenched case 

First we carry out a cluster expansion of the generating function (F x ) x . For 
this purpose we write the Hamiltonian H x ({S}) in the form 
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where 

Hi = RiSi', Ri = Hi + y ] r<fli> r<fii = ^ ] Xj a r (fj a ', (3-2) 

Uij = jtpiSi i(pjSj -\- KijSiSj . (3-3) 

We introduce here the parameters r (pi, which play a role of the effective 
field acting on the spin Si from the nearest neighbour at the site r. Summa- 
tion i n (|3-1| ) spans pairs of the nearest neighbour sites. We can write 
function F x in the form 

F x = Y^Fi + ]nQ . (3.4) 

i 

Here we use the notations 

Fi = In Zc, Zi = Sp 5i exp Hi ; 

Q = (expJ2Uij} P0 ; p = J\_pi\ Pi = Z^ 1 exp(flj) . (3.5) 

(ij) * 

The first term /Ci of a cluster expansion [ ||, ^] of In Q has the form 

In Q = ln(exp ^ Uij) Po /Ci = ^ In (exp L/y ) po = -z ^ F + ^ , 

(*i) * 

(3.6) 

where z = • 7Tjj is a number of the nearest neighbour sites and 
Fij = In Zij; Zij = Sp s . s . expHij; 

Hij = iKjSi + jKiSj + Kij SiSj] r R,i = Ri f(fi . (3-7) 

So the generating function (F x ) x restricted by the first term of the cluster 
expansion is obtained in the form 

(F x ) x = -z'J2{Fi)x + ^{Fi j } x (z' = z-l). (3.8) 

« (ij) 

Parameters r <Pia are found due to the minimization of the free energy 

d(F x ) x 



d r (fi 



. (3.9) 



It gives the following self-consistency equations for the average value of spin 

m£? = ((S ia ) H ) x 

= ((S ia ) Pi ) x , (3.10) 
((Sia) Pi ) X = ((S ia ) Pij ) x , (3.11) 



where j € 7Tj and 



Plj = ry • (3-12) 



Zi 3 
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Figure 1: Order parameters of pure and dilute (w\\ = cf) systems on plane 
square lattice (z = 4). 



m£? = (3.13) 



Relations ( 3,10| ), ( |3,11| ) can be written in the form: 

™V = <*£>> 
(4X = (F^) x - (3-14) 

Here the notations 

( F ial~ajx = JT. ( F i)x = (( S iai ' ' ' Sia n ) C pi )x 5 (3.15) 



/p( n ® d d d 

= ({Sicn ■ ■ ■ Si oln Sji3 1 ■ ■ ■ Sj^Yp.^x (3.16) 

are introduced. Below we give the explicit expressions for the quantities 
(Fi) x , {Fij)x an d their derivatives which will be used later: 

{Fi)x — ^ ' CgFiai Fia — hi 2 COsh K{ a , 

a 

I p( n ) \ — X X n p( n ) • 

\ r io.\ -a.nl x ~ Ua i a 2 ' ' ' u a n -xa n ^ai- r i ai ; 
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F ia = tanhre iQ ; 



= 1 - tanh 2 K ia ; 



,(3) 



_2F (1) F {2) ■ 

ia ia ' 



{Fij) x = WapFiajp; Fia,j/3 = In 2e Kal3 L ia j^ ; 

a/3 

£ia,i/3 = cosh(jK iQ + 1%^) + a a/3 cosh(jK iQ , - jK^); a a/3 = e~ 2K ^ • 

^iax-'-anJ/x — ai a 2 " ' ' °a„_ia n W aiP p ion,j/3 ' 

= (sinh(i«ia + + a Wj g sinh( J -« iQ - /L iadl3 ; 

p(20) 

iaj/3 

p(30) 
iaj/3 
,(n/) 



ia,j/3 J 



1 + a a/3 + 2a 0(/ 9siiih(2iKj ; a)) /L 



9p (10) p (20) . 
~* r ia,j/3 r ia,j/3 ' 



J ia,j[} 



1 - a 



a/3 . 



7 (H) 

ia,j/3 r2 

^iaj/3 



F (21) 

ia,j/3 



-2F {10) 

ia,j/3 ia,j/3 



(3.17) 



The expression ( 3.1l| ) contains iVzfi equations for the same number of vari- 
ables r (f>ia- At the case of unifo rm field («j Q — > K a ) parameters r (pi a lose 
their site dependence and ( 3.11 ) reduces to f2 equations for the same num- 
ber of fields and we obtain a well-known result of Bethe approximation and 
cluster variation method for average value of spin 



m 



(i) 



c a tanh K a = 22 



a/3 



P 



sinh(«4 + dp) + a a p sinh(K' Q - 
cos1i(k^ + R'a) + a a /3 cosh(K / Q — k'* 



(3.18) 



«a = K a + zip a ; K a = K a + z'ip a . (3.19) 

We can obtain CFs of any order with differentiation of ( 3.10| ), fl3,ll|) and it 
is an advantage of prese nted approach. 
It follows from (3.13) that 



_(2) _ 6™\a _ d ( F ia)x SKyy _ , (2) . (1) 



(3.20) 



Hereafter the summation over sort indices (Greek letters) is not written 
explicitly so the sum over 7 symbol is omitted in ( |3.20 ). We have just used 
the notation 



8 Hi, 



.i3 



(1) . 

J/3' 



_(1) _ 5 r ipi, 
riajP 



(3.21) 



(see a definition of Ki ( |3.2j )). It is convenient to write the relation ( p.20|) in 
the matrix form 

,(2)_ £.(2) fed) (3 _ 22) 

(3.23) 



where the matrix notation is introduced: 



a/3 



m- 



(2) 



ia,jf3 ' y i 



^(2) 



a/3 



(F {2) ) 



= (1) 



a/3 



,(1) 
v ia,i/3 
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In order to calculate the quantities r'-ffajp it is necessary to use the relation 
( 3.1 4| ) . Differentiating it with respect to Kjp we obtain: 



F (2) k (1) 



F, 



(20) g(l) 



(11) .&(!) 
rj 



where 



i(i) 



i(l) 



(20) v . 

ia/3,r/X' 



a/3 



ai3 



\ r ia,r/3/ x ' 



r ^ia,j/3 



(3.24) 

(3.25) 
(3.26) 



Taking into account (3.26), the equation ( |3.24| ) can be rewritten in the 
following form: 



p(20) -(1) , *=(1) _ / 

r ir r^ij < r ir tWrj ~ \ 



F : 



(20) p(2)\ g(l) 



,(H) -(1) 



l? y ' _|_ 77^ ' K - 

ij ' ir rj 



Replacing the indices i <-> r in ( |3.27 ) we get the second equation 



F (11) 



Vii ' r r% Wrj ~ r ri "'ij ' \ r ri r r J ™rj 



(3.27) 



(3.28) 



which forms together with ( |3.27 ) the set of equations for the quantities 



rftij , iVrj ■ Excluding iip)!j' from this set and expressing the k\j' via 
(see (|3.22|) ) we have 



,(i) 



(i) 



(2) 



™( 2 ) + n- m (2) 



(3.29) 



where 



Fj, 



#(n) / r F (20) V 1 F (11) -F. 

ir I n y ri i 
ri \ ir I ir i 



-1 



-1 



(3.30) 



Later o n we sum the relation ( |3.29j ) over r € TTi and take into account 
( 3.22| ), (|3,21 ). It results in the following Ornstein-Zernike-type equation for 

the pair CF m 



(2) 



rn. 



(2) 
rj 



where 



Un 



(3.31) 



(3.32) 



re^i 



In the case of uniform field the matrices F are independent of site indices, 

(2) 

the CF rh^j depends on the difference of the sites' 1, 2 coordinates and we 
can solve the equation ( 3.31| ) carrying out the Fourier transformation 

„ (2) ^ V f 



m 



X) 



rhP^ 1 (rj 



n) 



(27T) 



D 



dq e 



(3.33) 
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Figure 2: Phase diagrams for dilute system on plane square lattice (z = 4) 

at different " nonexchange" interaction Vn). Annealed: I - paraphase, II - 
spin ordered phase, III - segregated alloy. 
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Figure 3: MFA results for the phase diagram of the system. J a p = J2j Jia,j/3i 
hi = EjUuji - 2iii,j2 + Ii2,j2) (see the Appendix). 
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V e is a volume of the elementary cell, D being the dimension of the lattice 
(space). The solution has the form 



+ 



^(PWUfM) 1 - PW) 1 , (3.34) 



+ (tt(O) - Tr(g) ) lF ( > [ 



where 



vr(^) = eW'-«h vr(0)=z. 



For the hypercubic lattices one has 



D 



n(0) = 2^cos(gja) , 



i=i 



(3.35) 



(3.36) 



where a is the lattic e spa cing. 

The expression (3.34) in the limit c\ = 1 yields the TCA result for the 
ideal (one-sort) system [ [?]]. Taking q = in ( |3.34| ) one gets the earlier 
result [ ||, |6] . The known exact formula for t he one-dimensional system in 
zero external field [|l| also follows from ( 3.34| ) . 



(3) 

We will obtain the expression for the 3-site CF m\£ 2/ a 37 differentiating 
equation ( p. 31 ) which can be written in the form 



tt ( 2 ) x x n ( 2 ) 



(3.37) 



rG7ri 



(here free site indices are denoted with numbers and summation over re- 
peated sort indices is implied). The differentiation with respect to 
yields 



,( 3 ) 



,(2) 



,(3) 



(3.38) 

where the point denotes the derivatives of corresponding quantities: 



la 




16 



Gla 



r5,3"f 



3 7 

lot-r 7 rS 

5G la<r s _ \ g 




til 



(3) 

la, 2/3,37 



5m 



(2) 

la,2/9 



3 7 

! ' 7 2f3 

m (3) 

3 7 




(3.39) 
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In ( [3.39 ) the threepoles are introduced which allows us to write ( 3.3£j) in 
more convenient form: 

7 imf 2 ' t^n 1 




(3.40) 



It should be noted that equality of site indices in threepoles and matrices 
implies also the symmetry with respect to corresponding sort indices, i.e. 




18 IS- 




(3.41) 



3e 3e 

Differentiating expressions ( glSOD and (|3^32|) we find 




+ 



+ 




+ 



l l 




(3.42) 



When obtaining the formula ( 3.42 ) the identities of the following type are 
used: 



s If 



(2) 



8 k 



3" 



(if>) 



-1 



(5k 37 J 



(3.43) 



They can be obtained for any matrix A by means of differentiation of the 

identity AAr 1 = 1. 

Now w e sha ll obtain the expressions for the threepoles which have ap- 
peared in ( 3.42 ) 



"3o 



(3.44) 
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In the diagram form it can be represented as 




The other threepoles we find in the same way: 





1 r^l3 



r l« r 3 



where 




{(SiaSi^Sr^Jx , 
(Fla^lp^x = {{SlaS r j3S r y)p lr ) x ■ 



/ F (21) v 



(3.45) 



(3.46) 



(3.47) 



(3.48) 



The quantity which appears in (|3l6|) follows from ( gjgg) , ( gjp and 
can be expressed as 



±(1) _ ^(i) ~(i) _ p . (2) A (2) 
13 — 13 13 — ■*- J Xx''^X3 ^lv ''' , r3 



(3.49) 



Using (|Q5|)-( |3l6|) in ( gg| ) and taking into account (ggp , (|3T2p we get 

the expression for the threepole U involving only the quantities calculated 

(2) 

above, i.e. the matrices and threepoles F and the pair CF mly. The 

same procedur e is u sed to calculate the threepole G. Using the obtained 
expressions in ( |3.40 ) we get the equation for the 3-site CF: 




1 m 



(2) 
12 



+ 
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where 
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{g 1t 



V 




x i?(30) 

E 




E u i 



+ 



x i?(21) / \F( 21 ), 





t.i. 



■ E r 



+ 



(3.54) 



In ( 3.51| )-( p.54 ) "t.i." means twin to the previous (in parentheses) item, 

i.e. item where the matri x sub stitution E <-> G has been carried out. The 
last (in braces) item of ( 3.51| ) is an example of t.i. to the previous (in 
parentheses) one. 

Assuming uniform field and performing the Fourier transformation 





(fj - fi,f k - fi) 



(3.55) 



(27T) 



D 



7i e 



-7L^..-r,)_^_ [ d q 2 e -^(r k -n) \m( 3 )/ 



(27T) 




we can obtain from (3.50) the following expression for the 3-site CF 

m (2) (9i + 92) ( 




(91,92) 



V 




+ 



tt(9i) 



7r(9i + 92) 




5 7 +vr(gV X C 




+ 




\ m (2) (9i) 



(3.56) 



-1 



(3.57) 



/ m (2) (9"*2) 

where the equality which follows from ( |3 . 3 1| ) 

U - ir(qi + q 2 )G = (m^i^ + q 2 ] 
has been taken into account. 

4. Annealed case 



In this case we shall introduce into Hamiltonian additional variation pa- 
rameters ip. The quantity r tpi a plays the role of effective field acting on the 
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Figure 4: Temperature dependencies of spin-spin and sort-sort CFs' q = 
Fourier-transforms and of specyfic heat (in units of ks) of pure and dilute 
(ci = 0.6) systems on plane square lattice (z = 4). 
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spin of sort a in the site i from its neighbour (of any sort) in the site r. 
H({S,X}) takes the form 

H({S,X}) = J2Hi + J2 W V' 

* (v) 

where 

H{ — K{Si + [J,i, /ij — /ij -|- ^ r^i'i r^Pi — ^ Xi a r 1pi a 

Wij = - j<fiSi - iftjSj + KijSiSj - jipi - iipj + Vij . (4.2) 

Further, the technique of the previous section leads to the generating func- 
tion beeing (within TCA) of the following form: 

= + > ( 4 - 3 ) 

* (y) 



where 



^• = ln^; Ztj = Sp 5 . s .e^ ; (4.4) 



The conditions 



Hij — j^iSi + iKjSj + KijSiSj + 

+ jMi + iflj + ; (4.5) 
jfii = fii - jipi = ^ + ri>i- (4.6) 



dT dT 

0; = (4.7) 



lead to 



d rifia d r 1pi, 

(Sia) = {Sia) Pi = {Si a ) pi . , (4.8) 

(X ia ) = (X ia ) p% = (X ia ) pi . , (4.9) 

where 

Pi = Z~ l esxpHf, pij = Z^ 1 expiJ^-; j € m , (4.10) 
i.e. again equality of t he o ne-site and intracluster unary CFs is im plied . 



We note that equations (4.9) like the ones for the chemical potential (2.10) 
are not independent due to the identity 

= 1 • (4.11) 



Therefore the number of independent fields r ipi a and chemical potentials 
fiia decreases. We put 

r &n = 0; Hm = • (4.12) 
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In order to shorten following formulae let us intoduce a notation S 
{S,X}, i.e. 

Si a , if a < Q 



X ia _Q, if a > Q 



and n = {k,m}, (f = {if,i>}, 



M, 



(n) 



(Sj lQl • • • Si nCtrl 



With this new notations (^8) and (|4.9| ) can be written as 



F({k}) (4.13) 



M 



(i) 



Pij 



(4.14) 



The same technique as for quenched case yields a result for a pair CF 
in the same form: 



(m (2) (€>) 



where 

K') 



-i 



(l-z)(^ 2 )) 1 + z(^ (20) +^ (11) ) X + 



-1 



+ (vr(O) - Tr(q)) (f^ (f^) 1 ^ 2 °) - (445) 



a/3 5Ki a dK 



(SiQ,Sj^) c 



(SiSj) c (SiXjY 

(X-iSj) c (XiXj)' 



(4.16) 



_ (4-17) 
In equation ( 4. 15| ) the sort indices lay in the restricted interval a = 1 • • • 2f2 — 
1 due to ( 4, lip , ( [4.12 ); and the matrices ( |4.17| ) have the dimension (2f2 — 
1) x (2f2 — 1). The exact result [ |H f or the one-dimensional system (z = 2) 
is a partial case of the formula (4.17). 
The relation ( [4.14 ) is equiva 



ent to 

SPSiJCjPij = Pi 



and it yields 



(nO) 



■i J 



(n) 



(4.18) 
(4.19) 

(H) . 
(4.20) 

An expression for the 3-site CFs also remains the same form as in the 
quenched case 



Therefore Q4.15] ) can be written in more simple form 
M^(q) = (l - (l + z'f - *•($)/) ; / = (^ 2 )) _1 T' 
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^ (a) (a)/m<?(5) 




0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
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quenched — . 




MFA - 




annealed — - 
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Figure 5: <f-dependence of the diluted system CFs. c\ = 0.6, V a p = 0. 
T/T c = 0.5, 0.7, 0.9, 1.1, 1.3, 1.7. q = q x = q y . 



QUENCHED AND ANNEALED ISING SYSTEMS 



19 



?r(<fi 






7T(?1 + <?2) 

with ^4, B, C, D given by (]3751|)-( [3T5^ ) and intracluster CFs 




i/3 



d 3 T 



?7 




■3/3 



J a 



?7 

ia x yi/3 




»> 7 jP 




dKiadKipdn 

9 hZicyd fcj Kijry 



(4.21) 



(SiaSjpSj-r)^ ■ (4.22) 



.n 



5. Discussion 



The great fluctuations of quantities under discussion make the results of 
"effective field" theories (like TCA and MFA) worse. Thus for the present 
model a great difference in interactions (e.g., K a p = K5\ a (sort 1 is diluted 
by noninteracting impurties)) is expected to worse the quality of TCA re- 
sults whereas great z, when interaction K ia jp couples many sites, improves 
these results. Taking into account that TCA gives exact results for the 
one-dimensional system (z = 2), the diluted system on plane square lat- 
tice seems to be the most difficult test for TCA, which can discover all its 
shortcomings. That is why below we concentrate our attention at this case. 

Now let us consider a numerical investigation of the obtained results. In 
the quenched case we assume for simplicity w a p = c a cp (complete chaos: 
{Xi a Xjp) x = {Xi a ) x {Xjp) x ). The figures |2[ || show model's phase diagrams 
of two-sort system for different parameters of the Hamiltonian. Curie tem- 
perature T c and temperature of alloy spinodal segregation T s are found from 

(2) 

the condition M^(q = 0, T) — > oo. Apart from a complicated form of the 
TCA phase diagrams ^ with respect to those of MFA || one can also re- 
mark the rule, that great values of the quantities K = K\\ + K22 — 2ET12, 
J = Jn + J22 — 2 J\2 enhance both segregation and spin alignment, whereas 
great V\\ = V\\ + V22 — 2V12, hi = hi + I22 — 2/12 enhance segregation and 
therefore enhance spin alignment, if K > 0, J > 0, and supress alignment, 
if K < 0, J < 0. 
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Temperature dependence of order parameters m\ jc\ or m[ jc\ of pure 
and dilute systems are shown in the fig. [I]. TCA and MFA results are com- 
pared. The known exact result for the pure system on the plane square 
lattice is also depicted. For quenched system within TCA the quantity 

iji 1 ) = mi /cx stays less than 1 even at T = and vanishes at c\ < c p = — 
(below percolation point). This contradicts to MFA result = 1 at T = 
for all c\. Moreover, within MFA and = Jnm$ (q) / \c\T) depend 
on temperature, concentration and interaction strength only via T/T^ FA 
{T MFA = Cl j n / kB ), a (l) = a ^(T/T c MFA ), a® = a™ (q,T/T c MFA ). In 
TCA this takes place for the case uun = c\ only (all interacting spins con- 
stitute one infinite cluster). 

One can see in the figure |6| that within TCA pair CF of the quenched 

(2) 

diluted system (q) is not zero at T = and it yields infinite rise of the 

susceptibility at zero temperature x(o) — * 00 ixio) ~ m ii \Q)/T) whereas 
MFA gives x(o) ~> similarly to the case of the pure system. The TCA 
behavior of x(o) an d o^ 1 ) can be attributed to interacting spins out of the 
infinite cluster of interacting spins (finite size cluster effect [p|). Indeed, one 
can find that susceptibility of the isolated spins, which are always present 
in diluted system, diverges at T = 0. Finite size cluster effects are observed 
also in other approximations (e.g., x(0j of the quenched diluted system in 
the effective field approximation of Kaneyoshi et.al. [ [|). Figure y shows, 
that significant differencies between annealed and quenched CFs appear at 
low temperatures. 

The ref. [ 10 1, where the pair CF of the pure system in paraphase within 
TCA was obtained, reports that x(o) has a maximum as a function of T 
above T c for a fixed q greater than <fo- We do not find any maximum for 

x{q), but instead we notice a maximum of m^(q) for all ir(q) G 0, pr^-J 

in the pure and diluted systems. It should be noted that MFA does not find 

any maximum for m^(q). 
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Figure 6: CFs of the ideal (ci = 1) and diluted (ci = 0.6) systems. 



22 R. R. Levitskii, S. I. Sorokov, R. O. Sokolovskii 



References 



[1] B. Ja. Balagurov, V. G. Vaks, R. O. Zajtsev. Statistics of one- 
dimensional model of solid solution. / / Fiz. Tver. Tela, 1974, vol. 16, 
No 8, p. 2302-2309. (In Russian). 

[2] W. Y. Ching, D. L. Huber. Monte-Carlo studies of the critical behavior 
of site-diluted two-dimensional Ising models / /Phys. Rev. B, 1976, vol. 
13, No 7, p. 2962-2964. 

[3] E. Stoll, T. Schneider. Ising ferromagnet with quenched impurities //J. 
Phys. A, 1976, vol. 9, No 7, p. L67-L70. 

[4] V. G. Vaks, N. E. Zein. On the theory of phase transitions in solid 
solution. //J. Exp. Teor. Fiz., 1974, vol. 67, No 9, p. 1082-1100. (In 
Russian). 

[5] H. Sato, A. Arrott. Remarks on magnetically dilute systems. //J. Phys. 

Chem. Solids, 1959, vol. 10, p. 19-34. 
[6] R. R. Levitskii, S. I. Sorokov, R. O. Sokolovskii. Thermodynamics 

and relaxational dynamics of disordered Ising model. / / Lviv, 1993, 

43p. (Preprint / Ukr.Acad.Sci. Inst.Theor.Phys.; IFKS-93-17U). (In 

Ukrainian) . 

[7] I. R. Yukhnovskii, R. R. Levytskii, S. I. Sorokov. Thermodynamics and 
distribution functions of Ising model. Two-site cluster approximation. 
// Kiev, 1986, 33p. (Preprint / Acad. Sci. Ukr.SSR. Inst.Theor.Phys.; 
ITP-86-142R). (In Russian). 

[8] R. Kubo. Generalized cumulant expansion method. / /J. Phys. Soc. 
Japan, 1962, vol. 17, No 7, p. 1100-1120. 

[9] T.Kaneyoshi, I.Tamura and R.Honmura. Diluted Ising ferromagnet: Its 
physical properties. // Phys.Rev.B, 1984, vol. 29, p. 2769-2776. 
[10] S. Havlin, H. Sompolinsky. Cluster approximation for q-dependent cor- 
relations in magnetic and ferroelectric systems. / / Phys. Rev. B, 1982, 
vol. 25, No 9, p. 5828-5835. 

Appendix 

Here we calculate the thermodynamic potential and CFs of the system with 
long-range interaction. Let the Hamiltonian be 

L H({S}) = Y, ^iaSia + - Y, KiajpSiaSjfj + 
ia ia,j/3 

+ Yl U iaXia + 2 X/ ViaJflXiaXjp + A L H , (Al) 
ia iajfi 



jpSiaSjfi + - Y, liot.jpXiaXjp 



(A2) 



ia,j/3 



ia,j/3 



where the interactions J; 



ia,j/3, liajp are not restricted by the nearest neigh- 
bour sites. Consider first the quenched case. Within MFA the (A2) is taken 
in the following form 



A L H MFA = Y, Jia,jP L m^S jP - l -Y J -.i/3 Ma L mf p , (A3) 

ia,j/3 iajfl 
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where L m!^ = ((Si a )L H ) x is the one-site CF of the system, Jiajp 



Jij o 



The generating function L (F X ) X ({F}) of the 



system can be ex pres sed through the generating function (F x ) x ({n}) of the 
reference system (2.1), 



L (F x ) x (m) = (F x ) x ({k}) -W Jiajfi L m^ L m^ , (A4) 



2 ^ 



where 



hi — TjQ, -|- X{ a , 



ET 
•Jia,j/3 I'ljp ■ 



(A5) 



Here Aj a is the mean field caused by the long-range interaction. Taking into 
account the relation 



5k jp 5 



5T 



"' 30 jfl 



j/3 \ k-y 



the CFs of the system 



(A6) 



L (n) 
Tfl ■ ■ = • • • 

""iiai—i n a n rp. rp 



^({r}) 



are found to be 



L (i) (i) 



(A7) 

(A8) 
(A9) 



m 



,(2) 



,(3) 



— m ia,lS^ ls > n V rn nv,j0,k'y + 



(A10) 



where m.W, m,( 2 ), m( 3 ) are CFs of the reference system (2.1) in molecular 
field flABP . We note that in dA9|) , ( [A10|) the sum over repeated sort and site 
indices is implied. Fourier transformation ( |3.33| ), ( 3.55| ) allows us to write 
flA9p , flSISD in the form 



(All) 



L m (3) 



(9*1,92) 



f2^ (3) 

= (91 + 9*2)^(91 + 92)m^ 7 (gi + g 2 ) + 

+ m2 e (^l, 92) (^5/3 + ^(91) L "i (2) (9i)) x 

x (5, 7 + JeA(9"*2) L mf 7 ) (g 2 )) . (A12) 
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The solutions of equations ( All| ) , ( |A12[ ) read 



L rh^ 




(9i,92) 



(l-m( 2 )(g)J(g)) _1 




1 - J(q)m( 2 \q) 
(91,9*2) 



(A13) 

-1 



(l - J(q)rhW(q) 
In ( |A14[ ) the equality 

1 + J(q) L m^ (g) = (l _ j(q)mW (£)) _1 , 



-1 

(A14) 
(A15) 



which follows from ( |A13|) is used. 

In the case of annealed system the same scheme leads to similar relations. 



L T({T}) = T({k}) -W J ia>j p L M^ L M« 



ia,j/3 



L M {1) = M- (1) , 

ia tot , 



l m[ 1) 



(A16) 
(A17) 

(A18) 



where T = {T, u}, 

and sums over sort indices run within interval 1 • • • 20. Taking into account 
{Xin) = 1 — J2a=l (Xia), (|A16| ) can be rewritten in the form 



(A19) 



where 1 < a, (3 < 20—1 and Jj 



ia,j(3 



Jij 

/, 



is a matrix of (20— 1) x (20 — 1) 



size with (0-1) x (0-1) submatrix (jij) = Iia,jp-ha,jQ-Iin,jp+IiQ,jQ- 
For the two-sort system (0 = 2) we have 



JiXJX Jil,j2 

Ji2,jl Ji2,j2 

hl,j\ — 2iilJ2 + Ii2,j2 



"13 ' 

Making use of the identity (Xi a Xjo) c = — J2pZi(Xi a Xj/3) c one obtains the 
higher CFs in the form 



M(2) (q) = ( 1 - M< 2 > (g) J(g) ) (q) , 



(A20) 




(9*1,92) 



i-M( 2 )(9V(9) 




1-J®M( 2 )(0) 
(91,92) 



-1 



1- 3{q)M^{q) 



-1 



(A21) 
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